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MATHEMATICAL REPRESENTATION OF CURRENT DENSITY PROFILES FROM ION THRUSTERS 


by ThaineW. Reynolds 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 


A concise and convenient way of expressing ion beam 
spreading characteristics is presented. The two- 
parameter equation, j(z,0)/[ j(z,O)] = ^ ^)1 , 

is shown to be capable of representing experimental 
beam current density profiles from several bombard- 
ment ion thrusters. This is a far-field expression. A 
calculation of the development of the profiles downstream 
of an extended source indicates that ion beam profiles 
should be taken at least eight thruster radii away to be 
representative of the far-field behavior of the beam. 

Parametric plots (involving various combinations of 
n and \) of thrust factor, ion beam envelopes, and 
back- sputtering due to ion beam ero sion iUustrate the 
use of the ion beam equation in beam-interaction calcu- 
lations. Experimental ion beams illustrated had thrust 
factors varying from 0.95 to near 1.0 . 

Introduction 

A description of the ion beam density profile of elec- 
tron bombardment thrusters is needed for both mission 
analysis and spacecraft design. The angular spread of 
the ion beam as it leaves the thruster affects placement 
of possible spacecraft components, such as solar cells 
and thermal control surfaces, to avoid the problems of 
the ion beam striking sensitive surfaces and causing 
sputtering erosion. A loss in thrust is also involved 
when a significant fraction of the beam ions have trajec- 
tories away from the axis. 

This report presents a two -parameter equation 
which has been found useful in expressing anal 5 rtically the 
ion beam profile of a number of bombardment ion 
thrusters for which experimental profile data were 
available. The procedure for determining the values of 
the parameters for fitting the equation to experimental 
data is described. Finally, several calculations involv- 
ing interactions of the beam profile are presented to il- 
lustrate its utility . 


Profile Equation 


Form of Equation 


Ion beam density profiles from electron bombard- 
ment thrusters have been described as having a rela- 
tively flat central core and an exponential dropoff with 
increasing distance from the beam centerline. An 
equation which exhibits these characteristics is 

1(z,g) ^ Q-[\ (1 ~ cos 0)]^ / 

j(z,0) 


J(z,0) current density in radial plane at axial distance 
z and radial angle 9 (The current density is 
taken in the plane normal to the z -direction 
(radial plane) though the particle trajectories 
are at angles to this plane. ) 
j(z, 0) current density in radial plane at axial distance 
z on centerline 

9 half-angle from centerline to point of measure- 
ment 

n, \ parameters to be determined 
(All symbols are defined in appendix A. ) 

This equation is deduced, then, simply from obser- 
vation that experimentally determined profiles can be 
reasonably well described with the appropriate choice of 
the parameters n and A. Equation (1) is analogous to a 
similarly derived expression for density profiles of 
supersonic jets expanding into a vacuum. (2) 

Figure 1 shows the typical shapes of profiles ob- 
tained for n— 1 and 2 and for several values of A.. A 
higher X for any given n indicates a more collimated 
beam. Similarly, a higher n for any given X indicates 
a more collimated beam. 

Equation (1) is a far-field e3q)ression. It describes 
the expandii^ flow as if the ion trajectories were straight 
lines which appear to originate from a point source. The 
current density thus would vary along any trajectory in- 
versely as the square of the distance from the source. 
The flow is completely defined as a function of the angle 
9 . Equation (1) must necessarily become inapplicable 
close to the thruster. The distance from the thruster 
within which the equation is not useful is discussed in 
appendix B. 

The total current at any axial location z must be 
constant and equal to the emitted current. Thus, 


y* 2irrj(z,fl)dr = jjiRq 


Substitution of Eq. (1) into Eq. (2) and rearranging yield 

y ,3, 

I 2 


Letting 


Eq. (3) becomes 


where as shown in sketch (a), 
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( 5 ) 

Designating the bracketed term in Eq. (5) as e(n,X), 


The other profiles of Fig. 3 were taken at distance 
ratios z/Rq varying from 4 to 23. 3. 

It can be seen that Eq. (1) is capable of representing 
the experimental, ion beam, current density profile data 
reasonably well. The greatest differences between ex- 
periment and equation values occur at the low-density 
values where the accuracy of determination is least. 


1 

c(n,X) = - 
2 

Eq. ( 5 ) can be expressed 

j 




% 

c(n,X) — 
z^ 


( 6 ) 


( 7 ) 


Finally, Eqs. (7) and (1) combine to give the com- 
plete current density relation at any location, 

V . 

V V ^ 

I II 

The part I of Eq. ( 8 ) gives the variation along the 
centerline; part II gives the radial variation. Values 
of c(n,X) are plotted in Fig. 2. 


Comparison of Equation with Experimental Data 


Some of the information from Fig. 3 is summarized 
in table I. 

Determination of Parameter Values 

An experimental profile can be readily checked 
against Eq. (1) to determine the appropriate values of 
n and X. The value of X is determined, independently 
of n, from the angle at which j(z, 0 ^)/[ j(z, 0 >] = e“^ 
(i. e. , 0. 368), 


X = 


1 

1 - co s 0 0 


(9) 


Then, from any other experimental point, the value of n 
is determined. For instance, at the angle 0 q at 
which j(z, 0 p j)/(j(z, 0 )) = 0 . 1 , 


0.834 

ln[X(l - cos 0Q j)] 


( 10 ) 


The radial profile will thus fit at the three values of 
angle 0, 0, 00 , and 0 q 


The applicability of Eq. (1) to describing experimen- 
tal ion beam profile data is shown in Figs. 3(a) to (g) . 

In these figures the experimental profile data cited(^» 
are shown by the data symbols. The solid lines are the 
plots of Eq. (1) for the particular values of n and X 
cited on each curve. The value of z/Rq, the ratio of the 
distance downstream that the experimental profile was 
taken to the thruster radius, is also shown. 

Two sets of experimental data for a given thruster 
have ion beam density profiles at two or more axial dis- 
tances (Figs. 3(d) and (e)). In Fig. 3 (d), the data at 
axial distance ratios z/Rq of 14. 8 and 25. 2 fall on the 
same curve, showing that the assumption of linear tra- 
jectories between the two distances was applicable to 
this case, hi Ref. 1, from which the data of Fig. 3(e) 
were obtained, the thruster size is not given. The ac- 
tual linear distance at which the several profiles were 
taken is therefore shown in Fig. 3(e), rather than the 
distance ratios z/Rq which are cited in the other plots. 
The profiles at axial distances of 16. 5, 21.5, 28.4, and 
41 . 4 centimeters are closely approximated by the single 
profile equation with n = 1. 19 and X = 85. The data of 
the closest profile (z = 11.5 cm) may be too close to the 
thruster to expect far-field behavior (see appendix B). 

The furthest downstream profile (at z = 53. 4 cm) shows 
considerable asymmetry and is plotted separately in 
Fig, 3(e) along with the same equation plot as in the 
other part of this figure. The profiles of Fig. 3(e) also 
support the assumption of linear ion trajectories in the 
far field. Reasons for the asymmetry of the one profile 
were not discussed in Ref. 1 . 


Application to Certain Problems 
Thrust Factor 


The ratio of the axial component of thrust to the 
thrust that would be obtained if all trajectories were 
directed along the axis is referred to herein as the 
thrust factor Fj. Thus, 

/'OO . 

/ m(r)v cos 0 2m: dr 

F-r ( 11 ) 

m(r)v27rr dr 



The velocity v is assumed to be constant for any tra- 
jectory and to be directed in the 0 -direction. Equation 
( 11 ) can be expressed, then, as 




2 



(19) 


Numerical integration of Eq. (13) resulted in the values 
of thrust factor shown plotted in Fig. 4. 

Ion Beam Envelope 

The fraction of the total ion beam enclosed within a 
given half-angle 0^ expressed as 


S = S(0) sec 01 

where a is the angle between the surface normal and 
the incident ion direction. 

For the variation of fix with the position variables 
z and r , use was made of Eq. (8) . 


I(total) 


/ r=z tan 9x n 


' e"I^(l “ <20® ^)] 27rr dr 


e”[^(l " oos 0 )] ” 27rr dr 


(14) 


Substituting Eqs. (4) and (6) pto Eq. (14) yield 


Wj) 

I(total) 


= 2c(n,\) 


/ 


exp^ 


yi + 


h dy 


(15) 


It can be seen in sketch (a) that for the circular tar- 
get a = 0 , and that for the cylindrical tank wall 
a = (tt/ 2) - 0. Also the arrival rate on the cylindrical 
wall is related to the arrival rate on the target at the 
same location by 




Tank wall 


(Xi(Z 


,R) 


Target 


tan 0 (20) 


By converting to the dimensionless variable x, where 


X = 


L. 

R 


(21) 


Numerical integration of Eq. (15) for values of n = 1 and 
2 and for several X-values yielded the curves shown in 
Figs. 5(a) and (b), for the fraction of the total beam con- 
tained within the angle 0i^. 

Return of Sputtered Target Flux to Thruster 

The testing of ion thrusters in vacuum tanks causes 
sputtering of the facility walls where the ion beam im- 
pinges. The return of this sputtered material to the 
thruster face and condensation thereon was shown to af- 
fect the operation of glass-coated accelerator grids. 

In order to get some estimate of the amount of back- 
sputtered material to expect, the following model was 
used (sketch (b)), which made use of the proposed beam 
profile Eq. (1); 

The returning (back sputtered) flux rate ft 2 depends 
on the local emission rate of sputtered material and 
the view factor dFx^ 2 from the point of sputtering to 
the thruster center, (The variation across the thruster 
face will not be more than a few percent as long as the 
sputtered surface is about five or more thruster radii 
away. ) 

M2= / (19) 

^rea 

Assuming cosine -law distribution of the sputtered 
material and cylindrical symmetry, the back flux from 
the target and the wall can be written separately and 
summed. (View factors defined in appendix A are de- 
rivable from relations in Ref. 8 . ) 



The sputtered emission rate depends on the in- 
cident ion rate px and the sputtering ratio S, 

i^x =iiXxS (18) 

In turn, S depends on incident ion energy, target 
material, and angle of incidence. The angular variation 
was assumed to follow the relation(9) 


and using relations (4), (8), and (18) to (21), Eq. (17) 




Term I of Eq. (22) is the back flux from the circu- 
lar target. Term II is the back flux from the cylindrical 
wall. The two parts of Eq. (22) are plotted separately 
in Figs. 6 and 7 for values of n = 1 and 2 and for sev- 
eral values of 

In Figs. 6(a) and (b), for example, the ordinate pa- 
rameter 

V \RqJ 2S(0) 

is plotted. This parameter accounts for the flux ratio 
P 2 / from a circular target of radius r to the thruster 
exit center. 


The return flux from a target is obtained, then, as 




!^2\ 

^ /target 


R. 


V r2 2[S(0)1| 


2S(0) (23) 


target 


It can be seen, for instance, that for n = 1 and 
A - 50, no significant fraction of the total return flux 
from a circular target is obtained beyond about a 20^ 
half-angle. For n = 1 and A = 15, most of the return 


3 



flux is oonfiiUHi U> Mu' tarj^ol that is vvlthia a i>-valuo of \{0) rui*rt*iil tMuvU)stHl wilhin roue o{ hMlf-anf;it' 0, A 
alx>ut 45^'. For n - 2 and those saiuo A~va!iu>s, Iho 

cutoff angles are nearer ir>^^ and respeetively. l(total) total eurrtMit. A 


Figures 7(a) and (b) show Uu> return flux fn)m (lu> 
tank wall s to the tli ru ste r e xi t e e nte r . He rt^ t lu> o ixli na tt' 
gives the retux*n flux fi'om x - 0 out to any axial dis- 
tance x-z/li. (The reference radius H in the so plots 
is the tank radius, not the tlirusier radius. ) 'I'he rtdurn 
flux from the cylindrical wall is then 




2 2S(0) 
0 


wall 


Ho 

2S(l)) 




In this case it can be setni that lor n - 1 and \ = 50 
there is essentially no return flux from the wall between 
X = 0 and x = 1.5. The major jjortion of Uie return flux 
from the wall comes from the wall area between x 1.5 
and x = 10. For n- 1 and A = 15, return flux from 
the wall close to x ~ 0 is noted, but essentially no walN 
return flux comes from tlie wall beyond about x = S. 

For n = 2, the same A-values represent more colli- 
mated beams; hence, there is less waU-return flux fix)m 
similar facility dimensions. 


j(z,0) current density in nulial plane at axial ills- 
tance z and radial angle A/em*' 

j(/., r) current ii<'nsitv In radial plane at axial dis- 
tune e z and radial distanee r, A/em*' 

j(z,0) eurrent dimsity In r.adial plane at axial ills- 
tanee z on center I in«v. A/ein*" 

J(d»A*o) eurrent density in tiu' radial ()lane at Uiruster 
exit (z 0) and at radial distanee r^^, A/cm’*’ 

j(0, 0) current density In radial plane at thruster exit 
(z 0) at eenterline, A/env“ 

j average eurrent density of thruster (eq. (H2)), 

A/cm*^ 


m(r) in ass flux In radial plane at radius r, 

kg/(cm^)(see) 


Concluding llemarks 

A two-paramotcr equation is shown to be capable of 
representing experimental ion beam density profiles of 
bombardment ion thruster exhausts. The equation is a 
far-fleld relation (eq. (1)) and therefore will not be ap- 
plicable too close to tlie thruster exit. A calculation of 
the development of the profiles tbwn stream of an ex- 
tended source indicates that ion beam density pix> files 
should be taken at least eight thruster radii away to bo 
representative of the far-field behavior of the ion beam . 
The equation was shown to represent profiles taken us 
close as four thruster radii downstream. However, 
these may not be representative of the fur -field patterns. 


n parameter of Kq. (!) 

U radius of facility vacuum tank (sketch (a)), cm 

raiUus of Uirustor, cm 
r radius, mn 

r^j radiuH In Umistcr exit plane, cm 

S sputtering yield, atoms/incidiavt ion 

S(0) sputtering yield at aormal Inciilcncc, 

atoms/incident ion 


The particular advantage of being able to expriiss 
the ion beam profile in a concise analytical relation has 
been illustrated by parametric plots of thrust factor and 
ion beam envelope curves. Experimental beam pix>files 
illustrated had thrust factors varying from ulx>ut 0 . 1)5 to 
nearly 1.0. 

Use of the beam profile equation to ealculate back- 
sputtering from vacuum chamber walls during ion 
thruster testing was illustrated. 

Appendix A 

Symlx)ls 

c(n, A) parameter defined in Kq. (IJ) 

dA^ differential area of thruster, cm^ 

dFj 2 diiTerential view factor, fraidion of emitted 

flux from surface I that intiircopts surfaiui 2 

I current, A 


v velocity 

X nondimcnsional a.xial distance, z/R 

y nondlmensional radial distance, r/z (« tan 0) 

z axial distance variable, imi 

Zrp distance of target fmm thruster c.xlt plane, cm 

(V angle between surfucii iKxrmal and incidiait Ion 

direction 

/J angle in thruster cxli plam^ (skct<!fv (ir)) 

0 angle from centerline of ion boam 

0^, angle at which j(/‘, 0(.)/j(z, 0) equal s 0 . :u;h 

0 Q ( iui g !c a t w li i c h |(z ^ ^ j ) / j ( z , 0 ) ih|U : i1 s 0 . I 

0 I iH)n<v half-angle (inclosing given fraction of 

total ion l)cam (cqs. |M) and (15)) 



X parameter of Eq. (1) 

ju^ particle arrival rate at surface in target plane, 
atom s/(cm^) (sec) 

jLtg particle arrival rate at center of thruster exit 
plane, atoms/(cm^)|sec) 

If average emitted ion rate from thruster, 
ions/(cm^)(sec) 

particle emission rate from target surface, 
atoms/(cm^)(sec) 

Appendix B 

Development of Far -field Profile from Extended Source 


In order to get some indication of the closest ap- 
proach to an extended source (the thruster) that a far- 
field expression such as Eq. (1) could be useful, a cal- 
culation of the profiles as a function of distance from 
such a source was carried out. 


If it is assumed that nothing influences the ion tra- 
jectories once they are beyond the neutralization plane 
of the thruster, the development of a particular far-field 
profile, expressible in the form of Eq. (1), would only 
come about because the local emission sources (for in- 
stance, the individual grid holes) exhibit this same type 
of profile. The amoimt of emission current from the 
individual sources may vary. In this particular example 
a distribution of current density aero ss the radius of the 
thruster exit (z = 0) is assumed to be 


j(0»ro) 

j(O.O) 



(Bl) 


That is, individual emitting holes are not assumed, but 
rather a continuous emission source over the thruster 
face is assumed. The magnitude of the emission cur- 
rent varies as Eq. (Bl) and the distribution locally var- 
ies as Eq. (1). 


For the profile of Eq. (Bl) the current density at 
the cent^line j(0,0) is related to the average current 
density j by 

j(0,0)- 2,167 

where 

T_ I 


A schematic showing the geometric relations is 
shown in sketch (c). 


(B2) 

(B3) 


The flux to point z, r from the differential area of 
the thruster is 


= c(n,\) ^ ° (B4) 

KO.V ' '^2 


By using 


cos 0! 


+ r^ - 2r^r co s /3 + 


(B5) 


and 


7T ^0 


J (^ > ^q) "■ 2.16 j CO si 2 


(B6) 


Eq. (B4) can be written 


■1 r) _ 2.16 c(n,\) 

i w2 


f cos(l^\ 

I 


exp 


yr§ + r^ - 2r0r cos jS + z^ 


rodrod/3 (B7) 


Equation (B7) was integrated numerically . Results 
for n = 1 and X = 20 are shown plotted in Figs. 8 
and 9. Figure 8 shows the variation of j(z, r)/j with 
distance for several 0 -values. Values of j(z,r)/j ap- 
proach closely those given by the far-field expression 
within about six to eight thruster radii downstream. 

The same calculations are presented in a different man- 
ner in Fig. 9 . This figure shows constant current den- 
sity lines. The near-field data (dashed lines) blend into 
the far-field values (solid lines) by about six to eight 
thruster radii downstream. 
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TABLE I. - COMPARISON OF ION BEAM PROFILES WITH EQUATION (1) 


Thruster identification 

Axial distance 
of profile, 

a/Ro 



Diameter, 

cm 

Type 

Figure 

Refer- 

ence 

n 

■ 

30 

Glass grid 

3(a) 

3 

10 

1 

24.3 

0.954 

5 

Glass grid 

3{b) 

4 

23,3 

.92 

40.2 

.966 

5 

Glass grid 

3(c) 

(a) 

4 

1. 28 

18.2 

.951 

5 

Two grid 

3(d) 

5 

14.8 

.9 

63.3 

.980 





25. 2 

.9 

63.3 

.980 

(b) 

(b) 

3(e) 

1 

— 


85 

.99 

15 

Two -grid 

3(f) 

6 

4 


34.6 

®,976 






HQI 


955 

30 

Glass grid 

3(g) 

7 

4.8 



.960 


^Unpublished data by J. H. Molitpr and H. J. King of Hughes Research Laboratory, Malibu, Calif. 


^Not given in reference. 
®Case HI of ref. 6. 
**Case I of ref. 6. 
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Figured - Value of parameter c(n,M. 
(Eq. {6U 


Current density ratio, j(z, 0)/j(z, 0) 


CO 

o 

1 ^ 

s£> 





(a) Data of reference 3. Glass -grid, 30- (b) Data from reference 4. Glass-grid, 5- 

centimeter-diameter thruster; n = 1; centimeter-diameter thruster; n = 0.9% 

X = 24. 3; z/Rg = 10. X “ 40. % z/Rg » 23. 3; experimentai profiie 


taken 58.4 centimeters downstream. (Data 



Angle from centerline of ion beam, 0, deg 


(c) Unpublished data by J. H. Molitor and (d) Data from reference 5. Glass-grid, 5- 

H. J. King of Hughes Research Laboratory, centimeter-diameter thruster; n"0.9; 

Malibu, Calif. Glass-grid, 5-centimeter- X=63.3. 

diameter thruster; n = 1. ffi; X » 18. 2. 

Figure 3. - Comparison of profile equation (1) with data. Solid line is equation for n, X values cited. 
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(e-1) For various axial distance z. (e-2) For axial distance 

of 53. 4 centimeters, 

(e) Data of reference 1. Thruster type and diameter not specified* 



(f) Data of reference 6. Two-^rid, 15- (g) Unpublished data by R. Bechtel 

centimeter-diameter thruster; z/Rq- 40. of Lewis, Glass-grid, 30- 

centimeter -diameter thruster; 
n«0.8; A»42.9; Z/RgM-S; 
thrust factor, Fj = 0. 960. 


Pimiro ^ ~ rnrirhiHeri 



-6103 



Angle from centerlineof Ion beam, S, deg 
(Wn-Z 

Figure 5. - Fraction of total beam enclosed within given half-angle 
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Dimensionless radial distance, y « r/Z| h tan 6 


0 5 10 15 20 25 30 35 40 45 

Angle from centerline of ion beam, 0, deg 

(b) n -2. 

Figure 6, - Return flux from circular target. 
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